Developmental Changes in Genome Replication Progression in Pluripotent versus Differentiated Human Cells.

DNA replication is a fundamental process ensuring the maintenance of the genome each time cells divide. This is particularly relevant early in development when cells divide profusely, later giving rise to entire organs. Here, we analyze and compare the genome replication progression in human embryonic stem cells, induced pluripotent stem cells, and differentiated cells. Using single-cell microscopic approaches, we map the spatio-temporal genome replication as a function of chromatin marks/compaction level. Furthermore, we mapped the replication timing of subchromosomal tandem repeat regions and interspersed repeat sequence elements. Albeit the majority of these genomic repeats did not change their replication timing from pluripotent to differentiated cells, we found developmental changes in the replication timing of rDNA repeats. Comparing single-cell super-resolution microscopic data with data from genome-wide sequencing approaches showed comparable numbers of replicons and large overlap in origins numbers and genomic location among developmental states with a generally higher origin variability in pluripotent cells. Using ratiometric analysis of incorporated nucleotides normalized per replisome in single cells, we uncovered differences in fork speed throughout the S phase in pluripotent cells but not in somatic cells. Altogether, our data define similarities and differences on the replication program and characteristics in human cells at different developmental states.

Incomplete reprogramming of DNA replication timing in induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) are the foundation of cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which could affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing, a process linked to both genome regulation and genome stability, is efficiently reprogrammed to the embryonic state. To answer this, we compare genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicate their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibits delayed replication at heterochromatic regions containing genes downregulated in iPSCs with incompletely reprogrammed DNA methylation. DNA replication delays are not the result of gene expression or DNA methylation aberrations and persist after cells differentiate to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and influence the quality of iPSCs.

Genome replication in asynchronously growing microbial populations.

Biological cells replicate their genomes in a well-planned manner. The DNA replication program of an organism determines the timing at which different genomic regions are replicated, with fundamental consequences for cell homeostasis and genome stability. In a growing cell culture, genomic regions that are replicated early should be more abundant than regions that are replicated late. This abundance pattern can be experimentally measured using deep sequencing. However, a general quantitative theory linking this pattern to the replication program is still lacking. In this paper, we predict the abundance of DNA fragments in asynchronously growing cultures from any given stochastic model of the DNA replication program. As key examples, we present stochastic models of the DNA replication programs in budding yeast and Escherichia coli. In both cases, our model results are in excellent agreement with experimental data and permit to infer key information about the replication program. In particular, our method is able to infer the locations of known replication origins in budding yeast with high accuracy. These examples demonstrate that our method can provide insight into a broad range of organisms, from bacteria to eukaryotes.

Emergence of replication timing during early mammalian development.

DNA replication enables genetic inheritance across the kingdoms of life. Replication occurs with a defined temporal order known as the replication timing (RT) programme, leading to organization of the genome into early- or late-replicating regions. RT is cell-type specific, is tightly linked to the three-dimensional nuclear organization of the genome1,2 and is considered an epigenetic fingerprint3. In spite of its importance in maintaining the epigenome4, the developmental regulation of RT in mammals in vivo has not been explored. Here, using single-cell Repli-seq5, we generated genome-wide RT maps of mouse embryos from the zygote to the blastocyst stage. Our data show that RT is initially not well defined but becomes defined progressively from the 4-cell stage, coinciding with strengthening of the A and B compartments. We show that transcription contributes to the precision of the RT programme and that the difference in RT between the A and B compartments depends on RNA polymerase II at zygotic genome activation. Our data indicate that the establishment of nuclear organization precedes the acquisition of defined RT features and primes the partitioning of the genome into early- and late-replicating domains. Our work sheds light on the establishment of the epigenome at the beginning of mammalian development and reveals the organizing principles of genome organization.

RIF1 regulates early replication timing in murine B cells.

The mammalian DNA replication timing (RT) program is crucial for the proper functioning and integrity of the genome. The best-known mechanism for controlling RT is the suppression of late origins of replication in heterochromatin by RIF1. Here, we report that in antigen-activated, hypermutating murine B lymphocytes, RIF1 binds predominantly to early-replicating active chromatin and promotes early replication, but plays a minor role in regulating replication origin activity, gene expression and genome organization in B cells. Furthermore, we find that RIF1 functions in a complementary and non-epistatic manner with minichromosome maintenance (MCM) proteins to establish early RT signatures genome-wide and, specifically, to ensure the early replication of highly transcribed genes. These findings reveal additional layers of regulation within the B cell RT program, driven by the coordinated activity of RIF1 and MCM proteins.

Repli-seq Sample Preparation using Cell Sorting with Cell-Permeant Dyes.

Replication timing is significantly correlated with gene expression and chromatin organization, changes dynamically during cell differentiation, and is altered in diseased states. Genome-wide analysis of replication timing is performed in actively replicating cells by Repli-seq. Current methods for Repli-seq require cells to be fixed in large numbers. This is a barrier for sample types that are sensitive to fixation or are in very limited numbers. In this article, we outline optimized methods to process live cells and intact nuclei for Repli-seq. Our protocol enables the processing of a smaller number of cells per sample and reduces processing time and sample loss while obtaining high-quality data. Further, for samples that tend to form clumps and are difficult to dissociate into a single-cell suspension, we also outline methods for isolation, staining, and processing of nuclei for Repli-seq. The Repli-seq data obtained from live cells and intact nuclei are comparable to those obtained from the standard protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Live cell isolation and staining Alternate Protocol: Nuclei isolation and staining.

Insights into common fragile site instability: DNA replication challenges at DNA repeat sequences.

Common fragile sites (CFS) are specific genomic regions prone to chromosomal instability under conditions of DNA replication stress. CFSs manifest as breaks, gaps, and constrictions on metaphase chromosomes under mild replication stress. These replication-sensitive CFS regions are preferentially unstable during cancer development, as reflected by their association with copy number variants (CNVs) frequently arise in most tumor types. Over the years, it became clear that a combination of different characteristics underlies the enhanced sensitivity of CFSs to replication stress. As of today, there is a strong evidence that the core fragility regions along CFSs overlap with actively transcribed large genes with delayed replication timing upon replication stress. Recently, the mechanistic basis for CFS instability was further extended to regions which span topologically associated domain (TAD) boundaries, generating a fragility signature composed of replication, transcription and genome organization. The presence of difficult-to-replicate AT-rich repeats was one of the early features suggested to characterize a subgroup of CFSs. These long stretches of AT-dinucleotide have the potential to fold into stable secondary structures which may impede replication fork progression, leaving the region under-replicated. Here, we focus on the molecular mechanisms underlying repeat instability at CFSs and on the proteins involved in the resolution of secondary structure impediments arising along repetitive sequence elements which are essential for the maintenance of genome stability.

Three-dimensional nuclear organisation and the DNA replication timing program.

In eukaryotic cells, genome duplication is temporally organised according to a program referred to as the replication-timing (RT) program. The RT of individual genomic domains strikingly parallels the three-dimensional architecture of their chromatin contacts and subnuclear distribution. However, it is unclear whether this correspondence is coincidental or whether it indicates a causal and regulatory relationship. In either case, the nature of the molecular mechanisms ensuring this spatio-temporal coordination is still unknown. Here, we review recent evidence that begins to uncover the existence of a shared molecular machinery at the core of the spatio-temporal co-regulation of DNA replication and genome architecture. Finally, we discuss the outstanding, key question of the biological role of their coordination.

Somatic mutations inferred from RNA-seq data highlight the contribution of replication timing to mutation rate variation in a model plant.

Variation in the rates and characteristics of germline and somatic mutations across the genome of an organism is informative about DNA damage and repair processes and can also shed light on aspects of organism physiology and evolution. We adapted a recently developed method for inferring somatic mutations from bulk RNA-seq data and applied it to a large collection of Arabidopsis thaliana accessions. The wide range of genomic data types available for A. thaliana enabled us to investigate the relationships of multiple genomic features with the variation in the somatic mutation rate across the genome of this model plant. We observed that late replicated regions showed evidence of an elevated rate of somatic mutation compared to genomic regions that are replicated early. We identified transcriptional strand asymmetries, consistent with the effects of transcription-coupled damage and/or repair. We also observed a negative relationship between the inferred somatic mutation count and the H3K36me3 histone mark which is well documented in the literature of human systems. In addition, we were able to support previous reports of an inverse relationship between inferred somatic mutation count and guanine-cytosine content as well as a positive relationship between inferred somatic mutation count and DNA methylation for both cytosine and noncytosine mutations.

Mutational signatures association with replication timing in normal cells reveals similarities and differences with matched cancer tissues.

Mutational signatures' association with replication timing (RT) has been studied in cancer samples, but the RT distribution of somatic mutations in non-cancerous cells was only minimally explored. Here, we performed comprehensive analyses of mutational signatures in 2.9 million somatic mutations across multiple non-cancerous tissues, stratified by early and late RT regions. We found that many mutational processes are active mainly or solely in early RT, such as SBS16 in hepatocytes and SBS88 in the colon, or in late RT, such as SBS4 in lung and hepatocytes, and SBS18 across many tissues. The two ubiquitous signatures, SBS1 and SBS5, showed late and early bias, respectively, across multiple tissues and in mutations representing germ cells. We also performed a direct comparison with cancer samples in 4 matched tissue-cancer types. Unexpectedly, while for most signatures the RT bias was consistent in normal tissue and in cancer, we found that SBS1's late RT bias is lost in cancer.

Replication Timing of Gene Loci in Different Cell Cycle Phases.

Replication of distinct genomic loci occurs at different times during cell cycle. The replication timing correlates with chromatin status, three-dimensional folding, and transcriptional potential of the genes. In particular, active genes tend to replicate early in S phase, whereas inactive replicate late. In embryonic stem cells, some early replicating genes are not yet transcribed reflecting their potential to be transcribed upon differentiation. Here, I describe a method for evaluating the proportion of gene loci that is replicated in different phases of cell cycle thus reflecting the replication timing.

Imaging Method Using CRISPR/dCas9 and Engineered gRNA Scaffolds Can Perturb Replication Timing at the HSPA1 Locus.

Fluorescence microscopy imaging of specific chromosomal sites is essential for genome architecture research. To enable visualization of endogenous loci in mammalian cells, programmable DNA-binding proteins such as TAL effectors and CRISPR/dCas9 are commonly utilized. In addition, site-specific insertion of a TetO repeat array, coupled with TetR-enhanced green fluorescent protein fusion protein expression, can be used for labeling nonrepetitive endogenous loci. Here, we performed a comparison of several live-cell chromosome tagging methods, including their effect on subnuclear positioning, expression of adjacent genes, and DNA replication timing. Our results showed that the CRISPR-based imaging method can delay DNA replication timing and sister chromatid resolution at certain region. However, subnuclear localization of the labeled locus and gene expression from adjacent loci were unaffected by either TetO/TetR or CRISPR-based methods, suggesting that CRISPR-based imaging could be used for applications that do not require DNA replication analysis.

Replication timing and transcriptional control: beyond cause and effect - part IV.

Decades of work on the spatiotemporal organization of mammalian DNA replication timing (RT) continues to unveil novel correlations with aspects of transcription and chromatin organization but, until recently, mechanisms regulating RT and the biological significance of the RT program had been indistinct. We now know that the RT program is both influenced by and necessary to maintain chromatin structure, forming an epigenetic positive feedback loop. Moreover, the discovery of specific cis-acting elements regulating mammalian RT at both the domain and the whole-chromosome level has revealed multiple cell-type-specific and developmentally regulated mechanisms of RT control. We review recent evidence for diverse mechanisms employed by different cell types to regulate their RT programs and the biological significance of RT regulation during development.

The evolution of the human DNA replication timing program.

DNA is replicated according to a defined spatiotemporal program that is linked to both gene regulation and genome stability. The evolutionary forces that have shaped replication timing programs in eukaryotic species are largely unknown. Here, we studied the molecular causes and consequences of replication timing evolution across 94 humans, 95 chimpanzees, and 23 rhesus macaques. Replication timing differences recapitulated the species' phylogenetic tree, suggesting continuous evolution of the DNA replication timing program in primates. Hundreds of genomic regions had significant replication timing variation between humans and chimpanzees, of which 66 showed advances in replication origin firing in humans, while 57 were delayed. Genes overlapping these regions displayed correlated changes in expression levels and chromatin structure. Many human-chimpanzee variants also exhibited interindividual replication timing variation, pointing to ongoing evolution of replication timing at these loci. Association of replication timing variation with genetic variation revealed that DNA sequence evolution can explain replication timing variation between species. Taken together, DNA replication timing shows substantial and ongoing evolution in the human lineage that is driven by sequence alterations and could impact regulatory evolution at specific genomic sites.

Single-copy locus proteomics of early- and late-firing DNA replication origins identifies a role of Ask1/DASH complex in replication timing control.

The chromatin environment at origins of replication is thought to influence DNA replication initiation in eukaryotic genomes. However, it remains unclear how and which chromatin features control the firing of early-efficient (EE) or late-inefficient (LI) origins. Here, we use site-specific recombination and single-locus chromatin isolation to purify EE and LI replication origins in Saccharomyces cerevisiae. Using mass spectrometry, we define the protein composition of native chromatin regions surrounding the EE and LI replication start sites. In addition to known origin interactors, we find the microtubule-binding Ask1/DASH complex as an origin-regulating factor. Strikingly, tethering of Ask1 to individual origin sites advances replication timing (RT) of the targeted chromosomal domain. Targeted degradation of Ask1 globally changes RT of a subset of origins, which can be reproduced by inhibiting microtubule dynamics. Thus, our findings mechanistically connect RT and chromosomal organization via Ask1/DASH with the microtubule cytoskeleton.

scRepli-Seq: A Powerful Tool to Study Replication Timing and Genome Instability.

Advances in "omics" technology have made it possible to study a wide range of cellular phenomena at the single-cell level. Recently, we developed single-cell DNA replication sequencing (scRepli-seq) that measures replication timing (RT) by copy number differences between replicated and unreplicated genomic DNA in replicating single mammalian cells. This method has been used to reveal previously unrecognized static and dynamic natures of several hundred kilobases to a few megabases-scale chromosomal units called RT domains. Because RT domains are highly correlated to A/B compartments detected by Hi-C, scRepli-seq data can be used to predict the 3D organization of the genome in the nuclear space. scRepli-seq, which essentially measures the copy number, can also be applied to study genome instability.

Computational workflow for integrative analyses of DNA replication timing, epigenomic, and transcriptomic data.

Temporal profiling of DNA replication timing (RT) in combination with chromatin modifications, chromatin accessibility, and gene expression provides new insights into the causal relationships between chromatin and RT during cell cycle. Here, we describe a protocol for in-depth integrative computational analyses of Repli-seq, ATAC-seq, RNA-seq, and ChIP-seq or CUT&RUN data for multiple marks at various time points across cell cycle and changes in their interrelationships upon an experimental perturbation (e.g., knockdown or overexpression of a regulatory protein). For complete details on the use and execution of this protocol, please refer to Van Rechem et al. (2021).

dcHiC detects differential compartments across multiple Hi-C datasets.

The compartmental organization of mammalian genomes and its changes play important roles in distinct biological processes. Here, we introduce dcHiC, which utilizes a multivariate distance measure to identify significant changes in compartmentalization among multiple contact maps. Evaluating dcHiC on four collections of bulk and single-cell contact maps from in vitro mouse neural differentiation (n = 3), mouse hematopoiesis (n = 10), human LCLs (n = 20) and post-natal mouse brain development (n = 3 stages), we show its effectiveness and sensitivity in detecting biologically relevant changes, including those orthogonally validated. dcHiC reported regions with dynamically regulated genes associated with cell identity, along with correlated changes in chromatin states, subcompartments, replication timing and lamin association. With its efficient implementation, dcHiC enables high-resolution compartment analysis as well as standalone browser visualization, differential interaction identification and time-series clustering. dcHiC is an essential addition to the Hi-C analysis toolbox for the ever-growing number of bulk and single-cell contact maps. Available at: https://github.com/ay-lab/dcHiC .

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability.

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.